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The method of consistent potentials is used to explain how a minimally coupled (classical) scalar 
field can suppress Mixmaster oscillations in the approach to the singularity of generic cosmological 
spacetimes. 
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In their long-term study of the approach to the singularity in generic cosmological spacetimes , Belinskii, Kha- 
latnikov, and Lifshitz (BKL) concluded that generic cosmological spacetimes approach the singularity as a (different) 
Mixmaster universe p] at every spatial point. The Mixmaster universe is characterized by "oscillations" from one 
' Kasner solution ^ to another. When BKL considered the influence of a classical minimally coupled scalar field on 

■ the approach to the singularity, they found that it can suppress Mixmaster oscillations [||. While this result is now 
' well-known, the mechanism by which it happens is not widely appreciated. I wish to address this issue here. 

^-H The role of the scalar field may be easily clarified using the method of consistent potentials (MCP) originally due to 

1—1 ' Grubisic and Moncrief |^ . The MCP has been applied to a variety of cosmological spacetimes to explain the nature 
^ ^ , of the (numerically observed) approach to the singularity in spatially inhomogeneous cosmologie s . The MCP first 
assumes that the approach to the singularity is asymptotically velocity term dominated (AVTD) and then looks for 
\ a contradiction. If the model is AVTD, it approaches arbitrarily closely to a Kasner solution with a possibly different 
^ ■ Kasner solution at every spatial point. For any model, an asymptotic velocity term dominated (VTD) solution may 
I ' be found by neglecting all terms in Einstein's equations containing spatial derivatives and taking the limit as r — > oo. 
^ . (The MCP includes the implicit assumption that strong cosmic censorship holds for these cosmological models — 

■ i.e. there exists a foliation labeled by r such that some curvature invariant blows up as r — > oo.) If the model is actually 
OO . AVTD, then substitution of the VTD solution into the full Einstein equations will be consistent — asymptotically, all 

■ terms neglected in obtaining the VTD solution will be exponentially small. For convenience, the MCP will be applied 
to the Hamiltonian whose variation yields the relevant Einstein equations since exponentially small (growing) terms 
in the Hamiltonian will yield exponentially small (growing) terms in Einstein's equations upon variation. Any terms 
which cannot be made consistent with the VTD solution indicate, within the MCP, that the model is not AVTD. 
Here I shall consider essentially the same models as in but with the addition of a scalar field. I shall explore the 
influence of a massless, minimally coupled scalar field in some detail. The possibly interesting dynamics which may 



o 



o 

I . result from exponentially coupled scalar fields will be discussed elsewhere. For each cosmology, I shall use the MCP 

^L' to show how the scalar field yields an asymptotically velocity term dominated (AVTD) approach to the singularity. 
. . , Consider first the primary Mixmaster model — the spatially homogeneous vacuum Bianchi Type IX cosmology . 

^ ■ In the presence of a spatially homogeneous scalar field, (/'(t), Einstein's equations are obtained from the variation of 

[ the Hamiltonian constraint 7i = where 

2n = -pl + pl+pl+pl + v{n,f]±) + e''''v{cj,) . (i) 

In Eq. (H), f7 is 1/3 the logarithm of the spatial volume and P± the anisotropic shears with pQ and p± their conjugate 
momenta. The scalar field momentum, p^, plays a crucial role in the approach to the singularity. The spatial scalar 
curvature appears in V{fl, P±) =^g '^R for '^g and '^R respectively the determinant and scalar curvature of the spatial 
metric while V((/)) is the scalar field potential. The key to the infiuence of the scalar field is the additional kinetic 
term in Eq. (|l|). Any scalar field coupling which produces such a term will yield the same effect as long as V(0) does 
not contain terms exponential in (j). Let us consider 

H = HK+nv (2) 

where Tix contains all the momenta and 

First, we note that, in these variables, the (strong curvature) singularity in these models occurs as — > —00. Thus, 
unless V(0) contains terms which asymptotically depend exponentially on Q (through asymptotic linear dependence 
of (/) on il), e^^V((/)) ^ as ^ —00. For now, we shall ignore this term. 
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The MCP requires us to assume that 7i = T-Lk- Variation of this Hamiltonian yields equations with the solution 



(4) 



(5) 



where v± — p±/\pn\ and = p^/\pn\ with the momenta all constant. The Hamiltonian constraint 7i = becomes 
Hk = or 



(6) 



First consider the vacuum case — = 0. Then we can write v± in polar coordinates (with unit radius) in the 
anisotropy plane. The minisuperspace (MSS) potential V is dominated by the first three terms on the rhs of Eq. ^ 
so that 



2V « e'^^~'^'^+ + 64^+4/3+ +4\/3/3- _|_ g40+4/3+-4V3/3__ 



Substitution of (Q), with = cos 9, v- — sin 6, into ^ yields 

2V ~ g-4|r2|(l+2cose) g-4|0|(l-cose-y3sine) _^ g-4|n|(l-cose+y3sine) 



(7) 



(8) 



Except for (the set of measure zero) 9 — {0, 2tt/3, 47r/3}, any (generic) value of 9 will cause one of the terms on the rhs 
of (||) to grow. For example, 1 + 2 cos 6* < will cause the first term to grow. This condition arises for cos6' < — 1/2 
or -27r/3 <9 < 27r/3. 

With the addition of the scalar field, > 0, so that 



vi + v'i = l-vl<l. 



(9) 



No term in Eq. m) will grow if we can satisfy simulatneously with Eq. 



1 + 2w+ > 0, 
1-V+- Viv- > 0, 
1-V+ + VSv- > 0. 



(10) 



Since 



^ = 1 is no longer required, solution of Eqs. (llOh is possible if < 1/2 and < 1/12 which can occur 
if 2/3 < < 1. Since pn > decreases at each bounce |12[, any initial value of will eventually yield > 2/3. 

In the absence of the scalar field, the growing term in Eq^^ (^ causes a "bounce" (see Q) which changes the values 
of v± according to the prescription first given by BKL Q. If the new values of v± satisfy (|lo|), there will be no 
further bounces so that the solution will approach (^-(|). We see, then, that the main function of the scalar field 
is to weaken the restriction imposed by the Hamiltonian constraint on the gravitational "kinetic energy." In Fig. 1, 
trajectories with and without the scalar field are shown in the plane. In the vacuum case, away from bounces, 

the trajectory must fall on the "Kasner circle" defined by v'j^+v'^ = 1. For a small initial value of p^, the non-vacuum 
model closely tracks the other until the amplitude of grows sufficiently large. At this point, the trajectory deviates 
noticeably from the Kasner circle. Shortly thereafter, v± are able to satisfy ( [l0| ) so there are no more bounces and 
the values of v± no longer change. (This causes the trajectory in the w-t -plane to end.) 
It is possible to retain the Mixmaster (oscillatory) behavior if V{(j)) in M) is such that 



^60 



(11) 



where a, ( > 0. A potential of this type will, according to the MCP, yield an exponentially growing term in (|^) 
unless = 0. But — means that the usual Mixmaster oscillations will occur. Thus, given (O), there will be no 
way to avoid oscillations. Coupling between scalar field and gravitational degrees of freedom in e^^^ '^{4') could lead 
to very complicated behavior. BKL, in fact, by coupling an electromagnetic field to a Brans-Dicke-like scalar field 
so as to produce a potential term like (|l^) where a and b are functions of the electromagnetic field, claimed to have 
restored the oscillations suppressed initially by the scalar field alone Q . Potentials exponential in the scalar field can 
also arise in string theory 1 1^] . 

The most complicated models to which the MCP has been applied are vacuum cosmologies with a single spatial U{1) 
symmetry and spatial topology [p^p7|,p^. The degrees of freedom are {x, z. A, cj} with conjugate momenta 
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{Px,PztPAiP,i^}- The variables are functions of spatial coordinates u and v and time t. Here the Hamiltonian 
constraint {H = 0) is 

= Hk + Hv (12) 
where Ti.K contains only momenta and 



e^' - T^e-^' ( . (13) 

-e4^ + (l-a;2) e^' + {I + xf j 



Its variation yields Einstein's equations. A transformation A ^ A — 2t will restore the explicit time dependence in 
( |l^ ) found in [P"5|-^t|. To apply the MCP, the appropriate generalization of (^) is required. This is found by solving 
the velocity term dominated (VTD) equations obtained by neglecting all terms in Einstein's equations containing 
spatial derivatives. In the limit as r ^ co, wc find 

uj^ujQ^ p = -iv^, r = r°, A = Ao - uat, pa = (14) 

where Vz, v^p, xq, p^, wo, , Aq, and ua > are functions of u and v but independent of r. (The sign of ua is fixed 
to ensure collapse.) In the VTD limit, the Hamiltonian constraint ( p^ becomes 

- \pI + \p' + \pI = 0. (15) 

To obtain this expression, it was necessary to assume v^p > Q and z)z > so that the exponential terms containing 
pI. e"^^ and e^"^ in Tix would be exponentially small. However, v^p > causes the term ie^~^'^e"''VawVba; in ( p^ ) 
to grow unless 

lim (A - 2z - 4:(p),r = -wa + "^v^ + iv^ < 0. (16) 



This condition is obtained by noting that Vz > implies z— >— ooasr— s-ooso that from (|T^ behaves as e 



-2z 



As was shown in 1 16 Ij] , the VTD form of the Hamiltonian constraint (|15|) gives = 4i;^ + Avz, so that (^^ cannot 



be satisfied with v^p > and Vz > 0. On the other hand, if either of these is < 0, either r^e*'^ or p^e*^ in (|l^ ) will 
grow. This leads to the prediciton (observed numerically) that the approach to the singularity is oscillatory. 

Just as in the spatially homogeneous case, in the presence of a scalar field (j){u,v,T), the AVTD limit of the 
Hamiltonian constraint becomes 

0=-IpI + Ip^ + IpI + IpI (17) 

where p^ is the scalar field momentum. It is now possible to simultaneously satisfy > 0, > 0, and — wa + + 
Av^ < since v\ = Av^ + 4u^ + from (^7|) can be made arbitrarily large with a sufficiently strong scalar field. At 
any representative spatial point, after some number of oscillations, the momenta will move into the range needed to 
make the remaining evolution AVTD at that point. Presumably, as r ^ oo, the model will become AVTD almost 
everywhere. 

It is possible to display this effect in a computer simulation of the full Einstein equations for U{1) symmetric 
cosmologies with and without a scalar field. For convenience, p^ is chosen to be constant in space. In Fig. 2, the 
evolution of (p toward the singularity at three representative spatial points is shown for both cases. It is clear that 
the scalar field suppresses the oscillations. In Fig. 3, Lp{u = v,t) is shown for both cases. Note that the formation of 
ever smaller scale spatial structure is suppressed by the scalar field. Although these particular simulations can only 
be followed to r « 40, the MCP predicts that, for the scalar field model, at the spatial points of Figs. 2a and 2b, since 
Vip = (to machine accuracy) and Vip > respectively, no further bounces will occur. At the spatial point of Fig. 2c, 
at least one more bounce to change the sign of Vu, is expected. 
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In models with two commuting spatial Killin g fi elds, the Hamiltonian constraint enters in an interesting way. 
Consider the magnetic Gowdy model discussed in [l^Jl^. Its variables {P, Q, A} and conjugate momenta {ttp, ttq, tta} 
depend on spatial coordinate and time r. Here, Einstein's equations are obtained by variation of a Hamiltonian H 
which is not the Hamiltonian constraint Ti.. 

^ = [^P + ^-''"^^Q + e"'" (P'e +e'''Q,l )] + B\,eM+^y' (18) 
where B measures the strength of the magnetic field. The variation of (^8|) with respect to Tr\ yields 

A,. = - 4^ K + e-'^nl + e'^^ {P,l W^Q.l )] + B^e^'+^y^ (19) 

which happens to be a rewriting of the Hamiltonian constraint Ti. = 0. The AVTD solution obtained from the variation 
of (|l^) neglecting spatial derivatives is 

P = VT, Q^Qo, \ = -v^T (20) 



where v = lim 7rp/(27rA) depends on but not on r. As was described in [|T^, substitution of ( pO| ) into (18) shows 
that Vi = e^'^^ng/nx will grow if d < 0, V2 = e^'^^^^ Q,g /nx will grow if u > 1 while V3 = _B2g(A+T)/2 .^j^ gj.^.^^ -j- 
< V < 1. Thus there is no value of v consistent with the AVTD solution ( po|) . 

Now consider a scalar field with momentum n^. Eq. ( p^ for A,t gains an additional term and becomes 

so that the AVTD solution is (since the equations for P and Q remain the same) 

P^VT, Q^Qo, X = -v^t~vIt (22) 

where = tt^/{27t\). The magnetic potential will now grow if < + < 1. The additional scalar field kinetic 
energy means that we can have < w < 1 needed to keep Vi and V2 small with + > 1 as is needed to keep 
V3 small. Thus it is expected that, eventually, oscillations will cease as v and enter the required ranges at almost 
every spatial point. 

The results of numerical simulations of the magnetic Gowdy models (on rather than the solv- twisted torus [Q) 
with and without a scalar field are shown in Figures 4 and 5. In the former, the evolution toward the singularity of P 
at three typical spatial points is shown. For the scalar field model, w > at all the spatial points so that no further 
oscillations are expected. Once again, the scalar field is seen to suppress the bounces and the growth of small-scale 
spatial structure. The entire evolution P(9,t) is shown in Fig. 5. The larger the initial amplitude of tt^, the more 
quickly the bounces will be suppressed. 

As in the spatially homogeneous case, the inhomogeneous models should continue to oscillate for scalar field po- 
tentials with an exponential form as in(pT|) almost everywhere. Although the mechanism described here has been 
examined only in particular systems using arguments based on the MCP approximation (reinforced by numerical 
simulations of the full equations), we expect the behavior discussed here to be valid generically. The bottom line is 
that the scalar field kinetic energy relaxes the restrictions imposed by the Hamiltonian constraint on the gravitational 
kinetic energy. Since it is these restrictions which lead to the Mixmaster oscillations, relaxing them will allow an 
AVTD approach to the singularity. 
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FIGURE CAPTIONS 



Figure 1. Comparison of Bianchi IX MSS trajectories for models with (dashed line) and without (solid line) a 
scalar field. The Kasner circle is shown. Each saved data point of the scalar field model is shown with a dot. The 
code described in is used to evolve spatially homogeneous Bianchi IX models. The initial scalar field momentum 
is = .001. Both models start from the same initial data (a). At (b), it becomes clear that the scalar field model 
trajectory cannot reach the Kasner circle. The sequence of Kasner solutions (indicated by values of v±) terminates 
at (c) for the scalar field model. The vacuum model continues indefinitely although only the first few bounces are 
shown. The trajectory first reaches (c) for logj^Q = 3.45. 

Figure 2. Evolution of (p toward the singularity for U{1) symmetric vacuum (solid line) and scalar field (dashed 
line) models at representative spatial points. The scalar field is modeled hy — 5 independent of space and time. 
The simulation consists of 128'^ spatial grid points and is similar to those described in p^ . 

Figure 3. Evolution of ip toward the singularity for U{1) symmetric (a) vacuum and (b) scalar field models. A 
diagonal line u = v in the u-v plane is shown for the same simulation as in the previous figure with < u, w < 27r and 
2 < T < 39. 

Figure 4. Evolution of P toward the singularity for a magnetic Gowdy model (see jlj]) with (dashed line) and 
without (solid line) a scalar field. The initial scalar field momentum is = cos 9. There are 1024 spatial grid points 
in this simulation. 

Figure 5. P over the entire 9-t plane (a) without and (b) with a scalar field for the simulation of Figure 4 where 
O<0-eo<2TT for 6*0 = -tt/5 and < t < 61. 
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